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Abstract 

Assuming 2>-v mixing, massive Majorana neutrinos and neutrinoless double-beta 
((PP)ov-) decay generated only by the (V — A) charged current weak interaction via 
the exchange of the three Majorana neutrinos, we briefly review the predictions for the 
effective Majorana mass |<m>| in (/3/3)oi/-decay and reanalyse the physics potential 
of future (/?/3)oi/-decay experiments to provide information on the type of neutrino 
mass spectrum, the absolute scale of neutrino masses, and Majorana CP-violation in 
the lepton sector. Using as input the most recent experimental results on neutrino 
oscillation parameters and the prospective precision that can be achieved in future 
measurements of the latter, we perform a statistical analysis of a (/3/3)(v-decay half-life 
measurement taking into account experimental and theoretical errors, as well as the 
uncertainty implied by the imprecise knowledge of the corresponding nuclear matrix 
element (NME). We show, in particular, how the possibility to discriminate between the 
different types of neutrino mass spectra and the constraints on the absolute neutrino 
mass scale depend on the mean value and the experimental error of |<m>| and on the 
NME uncertainty. The constraints on Majorana CP-violation phases in the neutrino 
mixing matrix, which can be obtained from a measurement of | < m > | and additional 
data on the sum of neutrino masses, are also investigated in detail. We estimate the 
required experimental accuracies on both types of measurements, and the required 
precision in the NME permitting to address the issue of Majorana CP-violation in the 
lepton sector. 
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1 Introduction 



There has been remarkable progress in the studies of neutrino oscillations in the last 
several years. The experiments with solar, atmospheric, reactor and accelerator neutrinos 
[1-7] have provided compelling evidence for the existence of neutrino oscillations caused by 
nonzero neutrino masses and neutrino mixing 1 . The latest addition to these results are the 
Super-Kamiokande (SK) data on the L/i?-dependence of the (essentially multi-GeV) //-like 
atmospheric neutrino events [9], L and E being the distance traveled by neutrinos and the 
neutrino energy, and the new spectrum data of the KamLAND (KL) and K2K experiments 
[10, 11]. For the first time the data directly exhibit the effects of the oscillatory dependence 
on L/E and E of the probabilities of //-oscillations in vacuum [12]. As a result of these 
magnificent developments, the oscillations of solar u e , atmospheric and P^, accelerator 
(at L ~ 250 km) and reactor u e (at L ~ 180 km), driven by non-zero v- masses and //-mixing, 
can be considered as practically established. 

The evidences for //-oscillations obtained in the solar and atmospheric neutrino and KL 
and K2K experiments imply the existence of 3-v mixing in the weak charged-lepton current: 

3 

VXL = ^ Ulj VjL, l = e,n,T, (1) 
j=i 

where vil are the flavour neutrino fields, z/jL is the field of neutrino Vj having a mass rrij and 
U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [13], U = £/pmns- All 
existing //-oscillation data, except the data of the LSND experiment [8], can be described 
assuming 3-u mixing in vacuum; we will consider this possibility in what follows 2 . 

The PMNS matrix can be parametrized by 3 angles, and, depending on whether the 
massive neutrinos Vj are Dirac or Majorana particles, by 1 or 3 CP-violation (CPV) phases 
[17, 18]. In the standardly used parametrization (see, e.g. [19]), t/pMNS has the form: 



PMNS 



/ C12C13 S12C13 

-S12C23 - Ci 2 S 2 3Sl3e l<5 C12C23 - Sl2S23Sl3e i<5 

V S12S23 - c 12 c 23 s 13 e tS -c 12 s 23 - s 12 c 23 s 13 e iS c 23 c 13 e tS 




where c^- = cos^-, Sij = sin^-, the angles 9ij = [0, 7r/2], 5 = [0, 2tt] is the Dirac CPV phase 
and CK21, CC31 are two Majorana CPV phases [17,18]. One can identify the neutrino mass 
squared difference responsible for solar neutrino oscillations, Ar/ig, with Am^ = m\ — m\, 
Am| = AiTij! > 0. The neutrino mass squared difference driving the dominant z/ M — > u T 
(pn — > v T ) oscillations of atmospheric (p^) is then given by |Am^| = (Am^l = |Am| 2 | 3> 
Amjj. The corresponding solar and atmospheric neutrino mixing angles, 9 & and 8a, coincide 
with 9 12 and ^23, respectively. The angle 9± 3 is limited by the data from the CHOOZ and 
Palo Verde experiments [20]. 



1 Indications for (/-oscillations were reported also by the LSND collaboration [8]. 

2 The interpretation of LSND data in terms of — * v e oscillations requires (Ato 2 )lsnd — 1 eV 2 . The 
minimal 4-neutrino mixing scheme, which could incorporate the LSND indications for — > v e oscillations, is 
strongly disfavored by the data [14]. The z/-oscillation explanation of the LSND results is possible, assuming 
5- neutrino mixing [15]. The LSND results are being tested in the MiniBooNE experiment [16]. 
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Thus, the basic phenomenological parameters characterising the 3-z/ mixing are: i) the 
3 angles # i2 , 23 , 9i 3 , ii) depending on the nature of massive neutrinos 1 Dirac (5) or 1 
Dirac + 2 Majorana (5, a 2 i, «3i) CPV phases, and iii) the 3 neutrino masses, mi, m 2 , m 3 . 
Getting precise information about the //-mixing parameters is of fundamental importance 
for understanding the origin of neutrino mixing (see, e.g. [21]). 

The existing neutrino oscillation data allow us to determine AttIq, |Am|J, sin 2 9 Q and 
sin 2 29a with a relatively good precision and to obtain rather stringent limits on sin 2 # 13 
(see, e.g. [5,10,14,22,23]). The data imply that Am| = Am 2 21 ~ 8.0 x 10~ 5 eV 2 , \Am\\ ~ 
2.2 x 10" 3 eV 2 , sin 2 6 Q ~ 0.30, sin 2 26 A ~ 1 and sin 2 9 13 < 0.05. The sign of Am|, as is well 
known, cannot be determined from the present (SK atmospheric neutrino and K2K) data. 
In the case of 3-z/ mixing the two possibilities, Am 2 ^^ > or Am^^j < correspond to 
two different types of z/-mass spectrum: 

- with normal hierarchy (or ordering), m x < m 2 < m 3 , Am\ = Am^ > 0, and 

- with inverted hierarchy (ordering) 3 , m 3 < m x < m 2 , Am^ = Am 2 2 < 0. 

Depending on the sign of Am^, sgn(Am^), and the value of the lightest neutrino mass, 
min(mj), the z/-mass spectrum can be 

• Normal Hierarchical (NH): 

mi«m 2 « m 3 , with m 2 = y/Am 2 Q ~ 0.009 eV and m 3 = yf\Am\\ ~ 0.047 eV; 

• Inverted Hierarchical (IH): 

m 3 <C mi < m 2 , with mi j2 = a/| Am\\ ~ 0.047 eV; 

• Quasi- Degenerate (QD): 

mi = m 2 = m 3 , with m x = m 2 = m 3 = m , m 2 ^> [Am 2 ^, m ^ 0.10 eV. 

The precision on the mixing angles 821, O23, #13, and on Am 2 x and jAm^l, can be sig- 
nificantly improved in future z/-oscillation experiments (see, e.g. [26-30]). The sign of Am| x 
can be determined by studying oscillations of neutrinos and antineutrinos, say, — > z/ e and 
V(i — ¥ v e , in which matter effects are sufficiently large. This can be done in long-baseline 
i/-oscillation experiments running both with neutrino and antineutrino beams (see, e.g. [31]) 
or in the neutrino mode only [32,33]. If sin 2 26i 3 ^ 0.05 and sin 2 623 ^ 0.50, information on 
sgn(Am| 1 ) might be obtained in atmospheric neutrino experiments by investigating the ef- 
fects of the subdominant transitions z/^( e ) — > and v^{ e ) ~^ ^e(» of atmospheric neutrinos 
that traverse the Earth [34,35]. 

The neutrino oscillation experiments, however, cannot provide information on the abso- 
lute scale of neutrino masses (or on min(mj)) and thus on the possible hierarchical structure 
(NH, IH, QD, etc.) of the neutrino mass spectrum. The oscillations of flavour neutrinos, 
ui — > v\i and v\ — > vy, 1,1' = e,\i,r, are insensitive to the nature — Dirac or Majorana — of 
massive neutrinos Vj\ they are insensitive to the Majorana CPV phases a 2 i )3 i [17,36]. If Vj 
are Majorana fermions, getting experimental information about the Majorana CPV phases 

3 In the convention we use (called A), the neutrino masses are not ordered in magnitude according to their 
index number: Ato^ < corresponds to to 3 < mi < to-2. We can also always number the neutrinos with 
definite mass, in such a way that [24] ui\ < to-2 < TO3. In this convention (called B) we have in the case of 
the inverted hierarchy spectrum: Ato|, = Am| 2 , Am\ = Am^. Convention B is used, e.g. in [19,25]. 



3 



in C/pMNS would be a remarkably challenging problem 4 [24,39-41]. 

Establishing whether Vj are Dirac fermions possessing distinct antiparticles, or are Ma- 
jorana fermions, i.e. spin 1/2 particles that are identical with their antiparticles, is of fun- 
damental importance for understanding the underlying symmetries of particle interactions 
and the origin of //-masses. Let us recall that neutrinos Uj with definite mass will be Dirac 
fermions if particle interactions conserve some additive lepton number, e.g. the total lepton 
charge L = L e + + L T . If no lepton charge is conserved, the neutrinos Uj will be Majorana 
fermions (see, e.g. [42]). The observed patterns of ^-mixing and of |Am^| and Am^ can 
be related to Majorana Uj and the existence of an approximate symmetry corresponding to 
the conservation of the lepton charge V = L e — — L T [43,44]. The massive neutrinos 
are predicted to be of Majorana nature by the see-saw mechanism of neutrino mass genera- 
tion [45], which also provides an attractive explanation of the smallness of neutrino masses 
and, through the leptogenesis theory [38], of the observed baryon asymmetry of the Uni- 
verse. Determining the nature (Dirac or Majorana) of massive neutrinos Uj is one of the 
fundamental problems in the studies of neutrino mixing (see, e.g. [21]). 

If neutrinos Vj are Majorana fermions, processes in which the total lepton charge L is 
not conserved and changes by two units, such as K + — > 7r~ + /i + + fi + + (A, Z) — ► 
(A, Z + 2) + fi~ , etc., should exist 5 . The only feasible experiments that at present have 
the potential of establishing the Majorana nature of massive neutrinos are the experiments 
searching for the neutrinoless double beta ((/3/3)o y )-decay (A, Z) — > (A, Z + 2) + e~ + e~ (see, 
e.g. [42,47,48]). Under the assumption of o^-decay generated only by the (V-A) charged 
current weak interaction via the exchange of the three Majorana neutrinos Vj (nij ^ 1 eV), 
the dependence of the (/3/3) 0! ,-decay amplitude A(/3f3) 0v on the neutrino mass and mixing 
parameters factorizes in the effective Majorana mass <m> (see, e.g. [42,49]): 

A{(3t3) Qv ~<m> M , (3) 

where Ai is the corresponding nuclear matrix element (NME) and | < m > | is given by 

\<m>\ = \ mi \U cl \ 2 + m 2 \U e2 \ 2 e ia21 +m 3 \U e3 \ 2 e ia31 \ . (4) 

If CP-invariance holds 6 , one has [50] a 2 i = kn, a 3 \ = k'n, where k, k! = 0, 1, 2, and 

77 21 EE e iQ21 = ±1, 7/31 = e™ 31 = ±1 (5) 

represent the relative CP-parities of Majorana neutrinos v\ and v-i, and v\ and z/3, respec- 
tively. As eq. © indicates, the observation of (/3/3)o y -decay of a given nucleus, and the 
measurement of the corresponding half-life, would allow a determination of | < m > | only if 
the value of the relevant NME is known with a relatively small uncertainty. 

The experimental searches for (/3/?)o!,-decay have a long history (see, e.g. [47,49]). The 
best sensitivity was achieved in the Heidelberg-Moscow 76 Ge experiment [51]: 

|<m>| < (0.35 - 1.05) eV, at 90% C.L., (6) 

4 The phases 021,31 can significantly affect the predictions for the rates of (LFV) decays fi — > e + 7, 
t — > [i + 7, etc. in a large class of supersymmetric theories with see-saw mechanism of i/-mass generation 
(see, e.g. [37]). Majorana CPV phases might be at the origin of baryon asymmetry of the Universe [38]. 

5 The existing experimental constraints on the | AL\ = 2 processes have been discussed recently in, e.g. [46]. 

6 We assume that mj > and that the fields of the Majorana neutrinos Vj satisfy the Majorana condition: 
C{vj) T = Vj, j — 1, 2, 3, where C is the charge conjugation matrix. 
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where a factor of 3 uncertainty associated with the calculation of the relevant nuclear matrix 
element [49] is taken into account. A similar result has been obtained by the IGEX collabo- 
ration [52]: |<m>| < (0.33-1.35) eV (90% C.L.). A positive signal at > 3cr, corresponding 
to |<m>| = (0.1-0.9) eV at 99.73% C.L., is claimed to be observed in [53]. This result will 
be checked in the currently running and future (/3/?)oi/-decay experiments. Two experiments, 
NEM03 (with 100 Mo and 82 Se) [54] and CUORICINO (with 130 Te) [55], designed to reach a 
sensitivity of |<m>| ~(0.2-0.3) eV, are taking data. Their first results read (90% C.L.): 

|<m>| < (0.7-1.2) eV [54], |<m>| < (0.2-1.1) eV [55], (7) 

where the estimated uncertainties in the NME are accounted for. A number of projects 
aim to reach a sensitivity to |<m>| ~ (0.01-0.05) eV [48]: CUORE ( 130 Te), GERDA 
( 76 Ge), EXO ( 136 Xe), MAJORANA ( 76 Ge), MOON ( 100 Mo), XMASS ( 136 Xe), CANDLES 
( 48 Ca), SuperNEMO, etc. These experiments, in particular, can test the positive result 
claimed in [53] and probe the region of values of | < m > | predicted in the case of IH and QD 
spectra [25]. 

In the present article we reanalyze the potential contribution that the future planned 
(/3/3)o^-decay experiments can make to the studies of neutrino mixing. The observation of 
(/5/3)oj/-decay and the measurement of the corresponding half-life with a sufficient accuracy, 
would not only be a proof that the total lepton charge is not conserved in nature (see, 
e.g. [56]), but might provide also unique information i) on the type and possible hierarchical 
structure (NH, IH, QD, etc.) of the neutrino mass spectrum [19,25], ii) on the absolute scale 
of neutrino masses [58], and iii) on the Majorana CP- violation phases [24]. We consider 
3-z/ mixing, assume massive Majorana neutrinos and (/J/^o^-decay generated only by the 
(V — A) charged current weak interaction via the exchange of the three Majorana neutrinos. 
As input in the analysis we use the results of recent studies of the precision that can be 
achieved in the measurement of the solar neutrino and CHOOZ mixing angles 8u and 6*13, 
and of the neutrino mass squared differences Am^ and lAmJ-J, on which |<m>| depends. 
The uncertainty in the measured value of |<m>|, which is due to the imprecise knowledge 
of the relevant nuclear matrix elements, is also taken into account. All relevant errors are 
treated in a statistically self-consistent manner. 

Our work is a continuation of earlier studies 7 (see, e.g. [19,25,41,58-63]). It is stimulated 
by the remarkable progress recently made in the experimental studies of //-oscillations [3,5, 
10,11] and by the recent analyses [26-30,64,65] in which the prospects for high precision 
determination of sin 2 #12, sin 2 #13, Am^, and jAm^l in future //-oscillation experiments have 
been extensively investigated. As a result of these studies the experiments that can provide 
the most precise measurement of the //-oscillation parameters |<m>| depends on, have been 
identified and a rather thorough evaluation of the precision that can be achieved has been 
made. In view of these developments a re-examination of the physics potential of the future 
(/3/3)oi/-decay experiments is both necessary and timely. 

The outline of the paper is as follows. In Section |21 we briefly discuss the present status 
of the determination of, and the prospect for improvements of the precision on, the neutrino 
oscillation parameters relevant to the analysis of (/3/3)oi,-decay experiments. In Section El we 
review the predictions for the effective Majorana mass |<m>| as a function of the lightest 

7 For an extensive list of references see, e.g. [57]. 
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neutrino mass and the type of the neutrino mass spectrum, taking into account present and 
prospective uncertainties in the neutrino oscillation parameters. In Section 0] we present the 
results of a quantitative investigation of the potential of a future (/?/?) o^-decay experiment 
based on a x 2 -analysis. We show our results as a function of quantities such as the observed 
mean value of |<m>|, its experimental uncertainty, and the uncertainty in the NME. We 
evaluate the possibility to obtain information on the lightest neutrino mass, the type of 
v-mass spectrum and Majorana CPV phases from a (/3/5)o^-decay experiment. In the latter 
case we take into account the constraint on the sum of neutrino masses E, which could be 
provided by cosmological observations, and investigate in detail the accuracies on |<m>| 
and S, required in order to probe Majorana CP-violation in the lepton sector. Finally we 
conclude in Sectional 



2 The Neutrino Mixing Parameters and |<m>| 

One can express [66] the two larger neutrino masses in terms of the lightest one, 
min(mj) = mo = mum, and of Am Q and Am^ 8 . Within the convention we use, in both 
cases of v-m&ss spectrum with normal and inverted hierarchy, one has: Am = Am^ > 0. 
For normal hierarchy, min(mj) = mi, Am^ = Am^ > 0, ma = (mf + Am|)3, and m.3 = 
[m\ + Am})i If the spectrum is with inverted hierarchy, min(mj) = m 3 , Am\ = Am^ < 
and thus mi = (m| + |Am^| — Am^ = (m| + jAm^l) 5 , mi = (m 2 + |Am^J) 3 . For 
both types of mass ordering, the following relations hold: \U C \\ 2 = cos 2 6* (l — sin 2 #13), 
\U e 2\ 2 = sin 2 # Q (l — sin 2 13 ), and \U c3 \ 2 = sin 2 #13, 9 Q = 12 . Thus, in the case of interest 
the effective Majorana mass |<m>|, eq. (J3J), depends in general on: i) Am\ = AmL 32 i, ii) 
@o — #12 an d Affig = Ar/^j, iii) the lightest neutrino mass m , iv) the mixing angle #13, and 
v) the Majorana CPV phases 0:21,31 • 

The best fit value and the 95% C.L. allowed range of |Am^| found in a combined analysis 
of the atmospheric neutrino 9 and K2K data read [5, 14]: 

I Am! I = 2.2 x 1(T 3 eV 2 , 

(8) 

|Ami| = (1.7-2.9) x 1(T 3 eV 2 . 

Combined 2-v oscillation analyses of the solar neutrino and KL 766.3 tyr spectrum data 
show [10,22] that Am| and 6 Q lie in the low-LMA region: Am| = (7.9±g;|) x 10~ 5 eV 2 , 
tan 2 6 & = (0.40ig;g?). The high-LMA solution (see, e.g. [67]) is excluded at ~ 3.3a. Maximal 
solar neutrino mixing is ruled out at ~ 6a; at 95% C.L. one finds cos2# > 0.28 [22], which 
has important implications for |<m>| (see further). In the case of 3-u mixing, the v e 
and v e survival probabilities, relevant to the interpretation of the solar neutrino, KL and 
CHOOZ data, depend also on #13 [21,68]. A combined 3-z/ oscillation analysis of these data 
gives [14,22,23] 

sin 2 # 13 < 0.027 (0.047), at 95% (99.73%) C.L. (9) 

8 For a discussion of the relevant formalism see, e.g. [19,57]. Notice that in [19] mo was used in the case 
of QD spectrum to indicate Too = mi ~ TO2 ~ TO3. Here we extend this notation to indicate the smallest 
neutrino mass for each type of spectrum. 

9 The current atmospheric neutrino data are insensitive to #13 satisfying the CHOOZ limit [5]. 
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Furthermore, such an analysis shows [22] that for sin 2 13 ^ 0.02 the allowed ranges of Am^ 
and sin 2 21 do not differ substantially from those derived in the 2-u oscillation analyses, and 
that the best fit values are practically independent of sin 2 13 < 0.05. The best fit values 
and the allowed ranges at 95% C.L. read [14,22]: 

Am 2 = 8.0 x 10" 5 eV 2 , sin 2 21 = 0.31 , 

(10) 

Am 2 x = (7.3 - 8.5) x lO" 5 eV 2 , sin 2 12 = (0.26 - 0.36) . 

Existing data allow a determination of Airig, sin 2 9 & and |Am^| at 3a with an error of 
approximately 12%, 24%, and 50%, respectively. These parameters can (and very likely will) 
be measured with much higher accuracy in the future. The data from phase-Ill of the SNO 
experiment [3] 10 could lead to a reduction of the error in sin 2 9 12 to 21% [29,30]. If instead 
of 766.3 tyr one uses simulated 3 ktyr KamLAND data in the same global solar and reactor 
neutrino data analysis, the 3a errors in Am 21 and sin 2 9 12 diminish to 7% and 18% [30] . The 
most precise measurement of Am 21 could be achieved [29] using Super-Kamiokande doped 
with 0.1% of gadolinium for detection of reactor v e [64]: the SK detector gets the same flux of 
reactor v e as KamLAND and after 3 years of data-taking, Am 21 could be determined with an 
error of 3.5% at 3a [29]. A dedicated reactor v e experiment with a baseline L ~ 60 km, tuned 
to the minimum of the v e survival probability, could provide the most precise determination 
of sin 2 9\2 [30,65]: with statistics of ~ 60 GW ktyr and systematic error of 2% (5%), sin 2 9n 
could be measured with an error of 6% (9%) at 3a [30]. The inclusion of the uncertainty in 
013 (sin 2 13 <0.05) in the analyses increases the quoted errors by (1-3)% to approximately 
9% (12%) [30]. The highest precision in the determination of |Am|J = |Ato 31 | is expected to 
be achieved from the studies of z/^-oscillations in the T2K (SK) [69] experiment: if the true 
IAtoIjjJ = 2 x 10~ 3 eV 2 (and true sin 2 23 = 0.5), the 3a uncertainty in |Am^| is estimated to 
be reduced in this experiment to ~ 12% [26]. The error diminishes with increasing |Am 31 |. 

In what regards the CHOOZ angle 9± 3 , there are several proposals for reactor u e experi- 
ments with baseline L ~ (1-2) km [27], which could improve the current limit, sin 2 13 < 0.05, 
by a factor of (5-10): Double-CHOOZ (in France), Braidwood (in the USA), Daya-Bay 
(USA-China), KASKA (in Japan), etc. The reactor 0i 3 experiments can compete in sensi- 
tivity with accelerator experiments (T2K [69], NOz/A [70]) (see, e.g. [26]) and can be done 
on a relatively short (for experiments in this field) time scale. 

Information on the absolute scale of neutrino masses can be derived in 3 H /3-decay ex- 
periments [71-73] and from cosmological and astrophysical data. The most stringent upper 
bounds on the v e mass were obtained in the Troitzk [72] and Mainz [73] experiments: 

m Pe < 2.3 eV at 95% C.L.. (11) 

We have m Dp = m 12j3 in the case of the QD z/-mass spectrum. The KATRIN experiment [73] 
is planned to reach a sensitivity of m De ~ 0.20 eV, i.e. it will probe the region of the QD 
spectrum. The Cosmic Microwave Background (CMB) data of the WMAP experiment, 
combined with data from large scale structure surveys (2dFGRS, SDSS), lead to an upper 
limit on the sum of the neutrino masses [74]: 

J^rrij = E < (0.7-2.0) eV at 95% C.L., (12) 

j 

10 During this phase the neutral current rate will be measured in SNO with 3 He proportional counters. 
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where we have included a conservative estimate of the uncertainty in the upper limit (see, 
e.g. [75]). The WMAP and future PLANCK experiments can be sensitive to S = 0.4 eV. 
Data on weak lensing of galaxies by large scale structure, combined with data from the 
WMAP and PLANCK experiments, may allow S to be determined with an uncertainty of 
5 ~ (0.04-0.10) eV [76]. Similar sensitivities can be reached by analysing the distortions in 
the Cosmic Microwave Background due to gravitational lensing in a future high sensitivity 
experiment [77]. 

3 Predictions for the Effective Majorana Mass |<m>| 

Given Am?,, |Am^J, 9 & and sin 2 13 , the value of |<m>| depends strongly on the type 
of the neutrino mass spectrum (NH, IH, QD, etc.) and on the Majorana CPV phases of the 
PMNS matrix, 021,31 (see eq. (p£J|). In what follows we will summarise the current status of 
the predictions for | < m > | . 

Normal Hierarchical Spectrum. In this case ra\ <C m 2 <C ra^, m = m 1; and one 

has [19] 



where we have neglected mi in eq. (|14j). Although one neutrino, v\, effectively "decouples" 
and does not contribute to | < m > | , the value of | < m > | according to eq. (fT4^) still depends 
on the Majorana CPV phase difference = 031 — 0^21. This reflects the fact that in contrast 
to the case of massive Dirac neutrinos (or quarks), CP- violation can take place in the mixing 
of only two massive Majorana neutrinos [17]. 

Since at 95% (99.73%) C.L. we have [14,22] y / ZW| < 9.2 (9.4) x 10~ 3 eV, sin 2 < 
0.36 (0.40), y/Am\ < 5.4 (5.7) x 10~ 2 eV, sin 2 13 < 0.027 (0.046), and the largest neutrino 
mass enters into the expression for |<m>| multiplied by the factor sin 2 13 , the predicted 
value of I < m > | is typically ~ few x 10~ 3 eV: for sin 2 13 = 0.03, one finds | < m > | ^ 0.005 eV 
(using the data at 95% C.L.). Using the best fit values of the indicated parameters (see 
eqs. (JH1) and (|T0|) ) and the same value of sin 2 #13 = 0.03, we get |<m>| ^ 0.0042 eV. 

The minimal value of |<m>| in eq. (jl4j) is obtained when there is a maximal com- 
pensation between the "solar neutrino" term, v/ Am?, sin 2 9 & cos 2 #i 3 , and the "atmospheric 
neutrino" one, a/ Am\ sin 2 013. At 95% (99.73%) C.L. we have ^/Am 2 A sin 2 013 ^ 1.5 (2.7) x 
10 -3 eV, while the "solar neutrino" term takes values in the interval (2.1 — 3.2) x 10 -3 eV 
((1.9 — 3.6) x 10 -3 eV). Thus, at 95% C.L. the "solar neutrino" term is larger than the 
"atmospheric neutrino" one and | < m > | is bounded from below. However, this may not be 
true considering the current 99.73% C.L. intervals of allowed values of the relevant oscillation 
parameters. 



<m> 




(13) 



(14) 
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It follows from eq. (JT3jl and the allowed ranges of values of Am^, Ami, sin 2 , sin 2 i3 , 
as well as of the lightest neutrino mass mi and CPV phases «2i,3i, that in the case of 
spectrum with normal hierarchy there can be a complete cancellation between the three 
terms in eq. (jlHjl . and one can have [58] |<m>| = 0. 

Inverted Hierarchical Spectrum. For IH neutrino mass spectrum, Am\ < 0, m = 
m 3 , and m 3 < mi = m 2 = a/| Ami| = a/ Am^. Using eq. (J1J) we find [19,24]: 



(cos 2 + e ia21 sin 2 ) cos 2 9 13 Jml + | Am|| + m 3 sin 2 i3 e 



(15) 



m 2 



+ | Ami | cos2 013 (l - sin 2 20 sin 2 ^) 5 , (16) 



Ami | cos 2 0i3 (l - sin 2 20 sin 2 ^) 2 , (17) 

where we have neglected m 3 sin 2 0i 3 in eqs. ()16j1 and ()17jl . The term m 3 sin 2 0i 3 can always 
be neglected given the existing data: even in the case where the spectrum is with partial 
inverted hierarchy and m| ~ |Ami|, the minimum of the sum of the other two terms 
in |<m>| satisfies a/ m| + | Am 2 A \ cos 20 cos 2 i3 ^> m 3 sin 2 0i 3 , since the data on and 
0i 3 imply cos20 ^> tan 2 13 . Even though in eqs. (fTBj) and (|T7|) one of the three massive 
Majorana neutrinos "decouples", the value of |<m>| depends on the Majorana CP- violating 
phase ai2i- It follows from eq. (fTTj) that 



| Ami | cos26> cos 2 0i 3 < |<m>| < y| Ami| cos 2 9 13 . (18) 

The lower and upper limits correspond to the CP-conserving cases «2i = vr and «2i = 0. 
Most remarkably, since according to the solar neutrino and KamLAND data cos20 ~ 0.40 
and cos 26*0 ^ 0.28 at 95% C.L., we get a significant lower limit on |<m>| exceeding 
10~ 2 eV in this case [25,58]. Using, e.g. the best fit values of |Ami| and sin 2 we find: 
|<m>| ^ 0.02 eV. The maximal value of |<m>| is determined by |Ami| and, according to 
eqs. (jHJ) and can reach |<m>| ~ 0.055 eV. The indicated values of |<m>| are within 
the range of sensitivity of the next generation of (/3/3) ^-decay experiments. 

Let us note that if Ami < 0' Le - ^ the neutrino mass spectrum is with inverted hierarchy, 
an upper limit on E = (mi + ni2 + m 3 ) < 0.10 eV would imply m 3 ^ 0.02 eV and therefore 



rn 



<C | Ami | . m this case the spectrum would be of the IH type and eqs. (fTTjl and (JTHJ) 
would be valid. 

The expression for |<m>|, eq. (fTTj) . permits to relate the value of sin 2 a 2 i/2 to the 
experimentally measurable quantities [19,24] |<m>|, Ami an d sin 2 20 : 

A sufficiently accurate measurement of |<m>| and of |Ami| and Q , could allow us to get 
information about the value of «2i- If, e.g. the data show unambiguously that |<m>| < 
\J | Ami I cos 2 #i 3 , that would imply ct2i 7^ 0. If in addition the data show that |<m>| > 
\J | Ami I cos 26*0 cos 2 0i 3 , one could conclude that c*2i takes a CP-violating value. 





<m> 


2 > 


1 


Ami 


cos 


4 013 J 


sin 2 20 
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Three Quasi-Degenerate Neutrinos. In this case m = m\ = m 2 = m 3 , wig 3> 
|Am^|, m ^ 0.10 eV. Hence, |<m>| is essentially independent of Am\ and Am^, and 
the two possibilities, Am\ > and Am\ < 0, lead effectively to the same predictions 
for |<m>|. The mass mo coincides with the v e mass m ie measured in the 3 H /?-decay 
experiments: mo = m.p e . Thus, mo < 2.3 eV, or if we use a conservative cosmological upper 
limit [75], mo = S/3 < 0.7 eV. The QD t^-mass spectrum is realized for values of mo, that 
can be measured in the 3 H /3-decay experiment KATRIN [73]. The effective Majorana mass 
| < m > | is given by 

|<m>| = m o |(cos 2 # + sin 2 # e ia21 ) cos 2 #i 3 + e ia31 sin 2 # 13 1 , (20) 



= m |cos 2 # + sin 2 # e ia21 | =m yi-sin 2 2# sin 2 ^. (21) 

Similarly to the case of the IH spectrum, one has: 

m cos 200 ^ |<m>| < m . (22) 

For cos 2^0 ~ 0.40, favored by the data, one finds a non-trivial lower limit on |<m>|, 
|<m>| ^ 0.08 eV. For values of the parameters allowed at 95% C.L. one has |<m>| ?t 
0.06 eV. Using the conservative cosmological upper bound on E we get |<m>| < 0.70 eV. 
Also in this case one can obtain, in principle, direct information on one CPV phase from the 
measurement of | < m > | , mo and sin 2 29 Q : 

9 ct 2 i / |<m>| 2 \ 1 . , 

sin 2 — = 1 - l - =-!- — 5 . 23 

2 I m 2 J sin 2 2^ V ; 

The specific features of the predictions for |<m>| in the cases of the three types of neu- 
trino mass spectrum discussed above are evident in Fig. ^ where the dependence of | < m > | 
on m = min(mj) for sin 2 9 Q = 0.31 and sin 2 9 13 = 0.01 and 0.03 is shown. The figures are 
obtained by including a 2a uncertainty in the predicted value of |<m>|. The uncertainty 
in |<m>|, a(|<m>|), has been calculated by exploiting the explicit dependence of |<m>| 
on the oscillation parameters Am|, Am\, sin 2 ^ and sin 2 #13 and assuming the following 
la errors (achievable in the future) in the determination of the latter: a(Am 2 & ) = 2%, 
a(\Am 2 A \) = 6%, a(sin 2 12 ) = 4% and two values of a(sin 2 13 ) = 0.004 and 0.008. The cur- 
rent best fit values of Am 2 & and |Am^| have been used. The Majorana CP- violation phases 
«2i and 0:31 were varied over all possible values they can take 11 . For the NH and QD (and 
interpolating) spectra, the regions within the black lines of a given type (solid, short-dashed, 
long-dashed, dash-dotted) correspond to the four different sets of CP-conserving values of 
the two phases a 2 i and 031, and thus to the four possible combinations of the relative CP 
parities (7721, r] 31 ) of neutrinos z/ 1;2 and v\ 3 . (+1,+1) solid, (—1,-1) short-dashed, (+1,-1) 
long-dashed, and (—1,-1-1) dash-dotted lines. If the spectrum is IH, the contribution to 
|<m>| due to m 3 can be neglected and the predictions for |<m>| become practically inde- 
pendent of Q3i (7731). In this case the regions delimited by the black solid (dotted) lines cor- 
respond to 7/21 = +1 (7/21 = —1). In the case of CP- violation all colored regions are allowed. 

n It follows from eq. (0J) that |<m>| is symmetric under the transformations 021.31 — > 27r — 0121.31 • This 
implies that it is sufficient to consider values of 021,31 in the range [0,7r] to cover all possible physical 
configurations for | < m > | . 
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Figure 1: The predicted value of |<m>| (including a prospective 2a uncertainty) as a func- 
tion of min(mj) for sin 2 6 Q = 0.31 and sin 2 #13 = 0.01; 0.03, and two different assumptions on 
the error on sin 2 #13. For the NH and QD (and interpolating) spectra, the regions within the 
black lines of a given type (solid, short-dashed, long-dashed, dash-dotted) correspond to the 
four different sets of CP-conserving values of the two phases a 2 i and a 3 i, and thus to the four 
possible combinations of the relative CP parities (7721, V31) of neutrinos z/ 12 and (+1, +1) 
solid, (—1, —1) short-dashed, (+1, —1) long-dashed, and (—1, +1) dash-dotted lines. For the 
IH spectrum, the regions delimited by the black solid (dotted) lines correspond to ^21 = +1 
{i]2i = — 1), independently of 7/31. The regions shown in red/medium-gray correspond to 
violation of CP-symmetry (see text for further details). 



The regions shown in red/medium-gray are the so-called "just CP- violation" regions [19]: 
an experimental point in these regions would unambiguously signal CP-violation associated 
with Majorana neutrinos. In the regions shown in blue/dark-gray and in green/light-gray 
it is not possible to distinguish between CP-violation and CP-conservation, because of the 
uncertainty implied by the errors on the oscillation parameters. 

The impact the prospective errors on Am\, Am^ and sin 2 9 Q have on the predictions 
for |<m>| is, in general, very small. More specifically, in the case of QD spectrum, the 
contributions of a(Am\), a^Am^) and ^(sin 2 ^) in a{\<m>\) can be neglected and only 
cr(sin 2 # ) induces an uncertainty in |<m>| which can be as large as few % if «2i ~ tt- 
For the IH type of spectrum, both adAm 2 ^) and <7(sin 2 # 12 ) are relevant and contribute to 
an overall a{\<m>\) ~ several %. Also in this case the error on |<m>| due to cr(sin 2 # 12 ) 
increases as varies from to tt. For NH spectrum, the dominant source of error is 
a(sin 2 #13). It affects significantly the predicted value of |<m>|. This explains the different 
allowed ranges of values for |<m>| obtained for a(sin 2 #13) = 0.004 and 0.008. Notice that 
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the impact of the errors is larger the smaller |<m>| is, i.e. when (0J31 — c^i) approaches the 
value 71. 

If the spectrum is with normal hierarchy (Am\ > 0), |<m>| can lie anywhere between 
and the currently existing upper limits, eqs. © and (J7|). This conclusion does not change 
even under the most favorable conditions for the determination of |<m>|, namely, even 
when |Am^|, Amg, 8 & and #13 are known with negligible uncertainty. 

The "gap" between the predicted values of | < m > | in the cases of IH and NH spectra 
allows us, in principle, to distinguish between these two types of hierarchical spectra [25,58]. 
Establishing, for instance, that |<m>| 7^ but |<m>| < 10~ 2 eV would imply, within 
the 3-neutrino mixing scheme with Majorana neutrinos under discussion, that the neutrino 
mass spectrum is with normal hierarchy, i.e. Am\ > 0. Depending on the value of mi, the 
spectrum could be either normal hierarchical (NH) or with partial hierarchy [19]. Obviously, 
such a result would rule out both the IH and QD spectrum. 

If the results in [53] implying |<m>| = (0.1-0.9) eV are confirmed, this would mean, in 
particular, that the neutrino mass spectrum is of the QD type. In this case, however, the 
measurement of |<m>| cannot provide information on the sgn(Am^). 

It should be clear from the preceding discussion that, depending on the measured value 
of |<m>| 7^ 0, the (/J/^o^-decay experiments may or may not provide information on both 
the type of v mass spectrum (NH, IH, QD, etc.) and sgn(Am^). If |<m>| ~ few x 
10~ 3 eV < 10~ 2 eV, both the type of the spectrum and sgn(Am^) will be determined. For 
a/| Am 2 A \ cos 2^0 < |<m>| < a/ | Am\ | , it would be possible to conclude that sgn(Am^) < 
only if mo ^ 0.02 eV, i.e. m 2 , <C |Am^|. In a relatively narrow interval of values of 
mo ~ few x 10~ 2 eV, for which m 2 , ~ |Am^J, one can have both Am^ < and Am\ > 0. 
In the latter case the v mass spectrum is with partial hierarchy. If |<m>| ^ 0.10 eV, the 
v mass spectrum is QD and the measurement of |<m>| will provide no information on 
sgn(Am 2 v ). 

Finally, if neutrino oscillation experiments show that Am\ < and therefore the v mass 
spectrum is with inverted hierarchy, while in (/?/3)oi,-decay experiments only the upper limit 
|<m>| < y/\Arrr\\ cos 26 & cos 2 #13 is obtained, that would mean either that there is a new 
additional contribution to the (/3/3) 0v -decay amplitude which interferes destructively with 
that due to the light Majorana neutrino exchange, or that the massive neutrinos Uj are 
Dirac particles. Similar conclusion could be made if, e.g., the KATRIN experiment shows 
that mo ^ 0.2 eV and correspondingly the v mass spectrum is QD, while (/3/3)oi/-decay 
experiments demonstrate only that the upper limit |<m>| < mocos26 ) holds. 

4 Analysis of the Implications of a (/3/3)oz/-Decay Half- 
Life Measurement 

4.1 On the NME Uncertainties 

If the (/9/3)o^-decay of a given nucleus is observed, it will be possible to determine 
the value of |<m>| from the measurement of the associated half-life of the decay. This 
would require the knowledge of the nuclear matrix element of the process. At present there 
exist large uncertainties in the calculation of the (/3/3)oi,-decay nuclear matrix elements (see, 
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e.g. [49]). This is reflected, in particular, in the factor of ~ 3 uncertainty in the upper 
limit on |<m>|, which is extracted from the experimental lower limits on the (/3/3)o^-decay 
half-life of 76 Ge. 12 For some nuclei (such as 100 Mo, 130 Te, 136 Xe), the uncertainties can be 
even larger. Recently, encouraging results on the problem of calculating the nuclear matrix 
elements have been obtained in [78]. A discussion of the problems related to the calculation 
of the (/3/3)oi/-decay NME is outside the scope of the present work. We would like to only 
note here that the observation of a (/3/3) 0t ,-decay of one nucleus is likely to lead to searches 
and eventually to observation of the decay of other nuclei. It can be expected that such a 
progress will help, in particular, to solve the problem of the sufficiently precise calculation 
of the nuclear matrix elements for the (ft /3)q u - decay [79]. 



4.2 The Method of Analysis 

The experimental observable in (/3/3)o^-decay is the decay rate r o b s measured with 
an experimental accuracy <7(T obs ). The observed decay rate has to be compared with the 
theoretically predicted rate 

r th = G|A^| 2 (\<m>\(x)) 2 , (24) 

where G is a known phase space factor and Ai is the NME. In eq. (I24J) x = (x osc , x ^) are 
the parameters determining | < m > | , which we divide into parameters measured in oscillation 
experiments, whose values we are going to use as input in the analysis, 

Xosc = (#12, 013, lAmlJ, Att4) , (25) 

and parameters that are, in principle, accessible by (f3 (3) 0u - decay experiments, 

x pp = ( m o, sgn(Am 2 1 ), a 21 , a 31 ) . (26) 

To investigate the potential to get information on the parameters x^* from the result of a 
generic (/3/3)o-decay experiment, we convert the observed decay rate and the experimental 
error into an "observed effective Majorana mass" and its error by 



l< m >r-^-^_ ffw = ____ ff( r ob ,), (27) 

where [A^ol is some nominal (theoretically predicted) value of the NME. If |<m>| obs > nopp, 
a positive (/3/?) 0l ,-decay signal is observed at the no C.L. Otherwise only an upper bound on 
| < m > | is obtained. The quantity opp defined in eq. (}2*T|) is a measure for the "accuracy" of 
the experiment. Then we construct a \ 2 i n the following way: 



£ \<m>\(x) — \<m> 
X 2 (x° g u ,F)= min ± - 9 — - 9 . (2t 

te[l/VF,VF] ° 2 pp + e< 



obs 2 



12 For the uncertainty on the NME for the (/3/3)o ! /-decay half-life of 76 Ge commonly a factor of 10 is 
adopted. Since |<m>| depends on the square-root of the half- life, typical values for the current uncertainty 
on |<to>| are factors from 3 to 4. 
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The parameter £ takes into account the uncertainty on the NME, and it is defined by 



\M\ 



(29) 



\M 



where \M\ is the unknown true value of the NME and \A4o\ is the nominal value used in 
eq. (|2*7j) to obtain |<m>| obs . The theoretical error cr t h in eq. (J2%|) takes into account the 
uncertainty implied by the errors on the oscillation parameters; it is calculated by 



where the index % runs over the four oscillation parameters given in eq. (|25jl . Sx l osc is the 
uncertainty on the parameter x l osc , and we have used the fact that to very good approximation 
the errors on the oscillation parameters are uncorrelated (see, e.g. [14]). 

Assuming that the value of the NME is known within a factor F > 1 , for given parameters 
X/s/3 we minimise the right-hand side of eq. (|28|) with respect to £, allowing £ to vary within 
the interval [1/y/F, VP], A perfectly known NME corresponds to F — 1. Note that in this 
way we do not introduce a probability weight for the NME; all values between \A4q\/\/F and 
Vf\m o| are treated on an equal footing. This procedure is similar to the "flat priors" used 
in unitarity triangle fits of the CKM matrix in order to account for theoretical uncertainties, 
see e.g. [80]. We have adopted this method, since it is not possible to assign a well defined 
probability distribution to the parameter £, and therefore, specifying a range for £ without 
imposing any further weight seems to be the most reliable procedure. Since the choice of the 
NME uncertainty factor F is subject to some arbitrariness we shall show results for various 
values of F. 

To combine a measurement of the o^-decay rate with a constraint on the sum of the 
neutrino masses E (obtained, e.g. from cosmological/astrophysical observations), we gener- 
alise eq. (|28jl in a straightforward way. To take into account the correlations between | < m > \ 
and E induced by the uncertainties on Am^ and Am| 1; we use the following covariance ma- 
trix in the x 2 -analysis: 



where T\ = £|<m>|, T2 = E, and a\ v = opp and a^ xp = 0£ are the experimental errors on 
|<m>| obh and E, respectively. 

4.3 Constraining the Lightest Neutrino Mass 

We start the quantitative evaluation of the physics potential of a (/?/5)o y -decay observa- 
tion by discussing the information that can be obtained on the absolute value of the lightest 
neutrino mass ttiq. Given an experimental result on |<m>| from a (/3/3)oi/-decay experiment, 
one can infer an allowed range for mo for each type of neutrino mass ordering. The results of 
such an analysis are shown in Fig. |2j For given values of | < m > | obs and its experimental error 
app, we minimize the \ 2 °f ec L- p8~|) with respect to the phases 021 and a&i, and calculate 




(30) 




(31) 



14 



the allowed range for m at 2a by using the condition x 2 ( m o) < 4. In Fig. |2]we adopted the 
best fit values for Am^, | A-rrz^ | and sin 2 12 (see eqs. (JHJ) and (fTUj) ). and sin 2 9 13 = 0. We 
have verified that the results hardly change if the values of the oscillation parameters are 
varied within the present 3a ranges. The dashed lines in Fig. 121 correspond to the current 
uncertainties of Am^, jAm^l, sin 2 9\2 and sin 2 #13, while the solid lines are obtained using 
the following prospective la errors: o^Am 2 ^) = 2%, adAm^l) = 5%, a (sin 2 6 1 2 ) = 3% and 
o"(sin 2 #13) = 0.002. By comparing the dashed and solid lines in Fig. |21 one observes that 
improving the accuracy of the oscillation parameters has only a minor impact on the results: 
we find only small improvements of the constraints on mo, while the qualitative behavior 
is unchanged. Consider first the case of |<m>| obs = 0.2 eV, shown in the right column 
of Fig. El In this case a positive signal should be established with high confidence by the 
next generation of (/3/3)oj,-experiments. If the experimental error in |<m>| obs is sufficiently 
small (app ^ 0.06 eV for NME uncertainty factor F < 3), i) the NH and IH spectra will 
be excluded and hence, the neutrino mass spectrum will be proved to be QD, ii) itlq will be 
constrained to lie in a rather narrow interval of values limited from below by mo ^0.1 eV, 
and iii) no information on sgn(Am 2 1 ) will be obtained. The uncertainty in the NME directly 
translates into an uncertainty in ttiq. 

In the case of an "intermediate" value of |<m>| obs = 0.04 eV shown in the middle 
column of Fig. |2l a lower and an upper bound on mo can be established for Am| x > 
if C/3/3 ^ 0.017 eV: 0.01 eV ^ m ^ 0.1 eV. In the case of Am 2 x < only an upper 
bound will be obtained: m ^5 0.1 eV. This result can be easily understood from Fig. [TJ if 
I < m > I is sufficiently large and app is small enough, the branch corresponding to the normal 
hierarchical spectrum extending to m = can be excluded. 

Consider finally the left column of plots in Fig. El corresponding to a very small value 
of |<m>| obh = 4 x 10~ 3 eV. For estimated typical values of app of the next generation of 
(/3/3) y-decay experiments and the mean value of |<m>| considered, only an upper bound 
on I < m > I can be established. It is clear that in this case one gets also only an upper bound 
on mo- Moreover, from the panel corresponding to a known NME (F = 1) one observes that 
for (ambitious) experimental accuracies, i.e. for app ^ 7x 10~ 3 eV, the case of sgn(Am| 1 ) < 
(inverted mass hierarchy) can, in principle, be excluded. This is a consequence of the 
lower bound on |<m>| for inverted ordering, which follows from the fact that cos 2$i2 is 
significantly different from zero (see eq. (JTHj) and the related discussion). However, if we 
take into account a possible uncertainty in the NME, the requirements on the experimental 
accuracy of |<m>| become exceedingly demanding (app ^4x 10~ 3 eV for F — 2), which 
renders the exclusion of the neutrino mass spectrum with inverted hierarchy remarkably 
challenging. Reducing the error to app = 10~ 3 eV would allow, e.g. for F < 2, to conclude 
that m ^ 0.02 eV and the neutrino mass spectrum is normal hierarchical. Establishing in 
an independent experiment that sgn(Am 2 1 ) < would imply in the case under consideration 
that there are additional mechanism(s) of (/?/3)o!,-decay [81] whose contribution to the (l3f3)o v - 
decay amplitude compensates partially the one due to the Majorana neutrino exchange. 

4.4 Determining the Type of Neutrino Mass Spectrum 

As is clear from the previous discussions, (/5/5)ov-decay experiments provide a unique 
possibility to obtain information on the type of neutrino mass spectrum, i.e. to distinguish 
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\<m>\ = 0.004 eV \<m>\ ' = 0.04 eV \<m>\ = 0.2 eV 




lo experimental accuracy on \<m>\ [eV] 
K\\\1 normal ordering V//A inverted ordering 

dashed: o(sin 2 9 13 ) = 0.016, a(sin 2 9 12 ) = 7.5%, o(Arn zl ) = 4%, a(Arn 31 ) = 13% sin 2 9 13 = 

solid: a(sin 2 9 13 ) = 0.002, a(sin 2 9 12 ) = 3.0%, a(Am 2I ) = 2%, o(Am\ x ) = 5% sin 2 9 12 = 0.31 

Figure 2: The reconstructed range for the lightest neutrino mass m at 2cr C.L. for normal 
(Amjj > 0) and inverted (Am^ < 0) mass ordering as a function of the la experimental 
error on |<m>| ob \ The results are shown for three representative values |<m>| obs = 
0.004, 0.04, 0.2 eV (columns of panels), and for fixed NME (first row), and an uncertainty of 
a factor of F = 2 and F = 4 in the NME (second and third rows). The figure is obtained 
using the current best fit values of Am^, lAmU and sin 2 12 (eqs. (jSJ) and (fTOjl). and 
sin 2 9\3 = 0. The dashed (solid) lines correspond to the present (prospective) uncertainties 
on the oscillation parameters. To the left of the dotted lines, a positive signal is obtained at 
2a, whereas to the right only an upper bound can be stated. 



between the NH, IH and QD spectra. As we have commented earlier, getting information 
on the possible hierarchical structure of the neutrino mass spectrum and on sgn(Am| 1 ) 
are different, although not totally unrelated, problems. In this subsection we elaborate 
further on the issue, since a determination of the neutrino mass spectrum is fundamental 
for our understanding of neutrino mixing. We investigate what conclusions can be drawn 
at the 2a C.L. on the neutrino mass spectrum from a result of a (/3/3) ^-decay experiment, 
characterised by the observed value of | < m > | and its experimental error. For given values 
of |<m>| ob& , a/3/3, the uncertainty F in the NME and a fixed neutrino mass ordering, we 
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observed l<m>l [eV] 




I I No information on the mass ordering Y///A Either IH or QD spectrum 

I I Inverted ordering excluded at 2a I I QD with no information on ordering 

To the right a signal is observed at 2a To the right the NH spectrum is excluded 

sin 2 9 13 = 0.03 ± 0.006, sin 2 12 = 0.31 ± 3%, Am^ = 8xl0~ 5 ± 2%, lAm^l = 2.2xl0" 3 ± 3% 



Figure 3: Information on the type of neutrino mass spectrum, inferred from data on (/3(3)q v - 
decay as a function of the observed |<m>| and its experimental error for three different 
assumptions on the NME uncertainty factor F (see text for details). 



minimize the x 2 -f un ction of eq. (|2*%|) with respect to mo,a2i and a^i- If the x 2 -minimum 
is smaller than 4, we conclude that this type of ordering is allowed. In addition we test if 
the "data" are consistent with negligible m , which implies a hierarchical spectrum (more 
precisely, we test whether x 2 ( m o = 0) < 4). The results of our analysis are shown graphically 
in Fig. El For values of |<m>| obs and &pp forming the dark shaded and white areas in 
the three panels, no information on sgn(Am| 1 ) can be obtained. The light shaded regions 
correspond to the case where sgn(Am| 1 ) < (inverted mass ordering) can be excluded. In 
agreement with the results presented in the previous subsection, we find that this is only 
possible for |<m>| obs < 0.01 eV and an experimental error well below 0.01 eV. To the right 
of the solid curve, the spectrum cannot be hierarchical for sg^Am 2 ^) > 0, i.e. the possibility 
mi <C m 2 <C m 3 is ruled out (at 2a). In the hatched region in this domain the "data" are still 
consistent with mo = for the inverted ordering, i.e. with an IH spectrum. Hence, if a result 
within the hatched region is obtained, we can conclude that either sgn(Am| 1 ) > and the 
spectrum is with partial hierarchy or of the QD type, or sgn(Am| 1 ) < and the spectrum 
is IH (m < 0.02 eV), or QD (m ^ 0.1 eV), or with partial hierarchy (m 2 , ~ |Am^|). This 
situation corresponds, e.g. to the panels of the middle column in Fig. or to the case when 
the lower branch at mo ^ 0.01 eV in the case of normal ordering (see Fig.QJ can be excluded. 
Finally, for sufficiently large values of |<m>| obs , corresponding to the white regions in Fig. El 
the spectrum is of QD type and no information on sg^Am 2 ^) can be obtained. 

Let us add that these results are stable with respect to variations of the oscillation pa- 
rameters within the present allowed ranges. In particular, they practically do not depend 
on the value of sin 2 #13: no significant changes appear if, instead of the rather large value 
sin 2 #i3 = 0.03 adopted in Fig. El smaller values are used. As discussed in Sectional cosmol- 
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Uncertainty in l<m>l from NME: 




observed \<m>\ [eV] 

Ix^NNNI Normal ordering Normal hierarchical I I \<m>\ and E inconsistent at 2a 

t///A Inverted ordering Inverted hierarchical 

sin 2 9 13 = ± 0.002, sin 2 12 = 0.31 ± 3%, Am^ = 8xl0" 5 ± 2%, lAm^l = 2.2xl0" 3 ± 3% 

Figure 4: Information on the neutrino mass spectrum from a combination of ((3/3)o u -decay 
data on |<m>| and cosmological data on E. These results are obtained for different as- 
sumptions about the errors in the determination of | < m > | and E (rows of panels) and the 
NME uncertainty factor F (columns of panels). We test whether the (/3/?)o;,-decay and cos- 
mological data are consistent with each other, with normal or inverted mass ordering, with 
normal ordering and m = (NH u-mass spectrum), or with inverted ordering and m = 
(IH spectrum), at 2a CL. The regions to the right of (above) the vertical (horizontal) dotted 
lines correspond to non-zero observed |<m>| obh (E obs ) at 2a. 

ogy provides a sensitive tool to constrain the sum of the neutrino masses E = mi + m 2 +m 3 . 
In the following we investigate what can be learned from data on (/5/5)o^-decay, combined 
with information on E from cosmology. We include the latter by assuming an "observed" 
value of the sum of the neutrino masses E obs with an experimental accuracy a-%- Obviously, 
if E obs > na-£,, cosmological observations would provide positive evidence for a non-zero E obs 
at the na C.L.; otherwise an upper bound is obtained. In Fig. 01 we show the results of a 
combined analysis of a (/3/3)o y -decay result with information from cosmology as a function of 
the observed values of | < m > | and E for two sets of representative experimental errors (up- 
per and lower rows) and different assumptions on the uncertainty from the NME. For given 
values of |<m>| obs and E obs we look for the x 2 -minimum for each of the two possible mass 
orderings. If a minimum is less than 4, we conclude that the corresponding mass ordering is 
consistent with the data at 2a. These regions are indicated by the hatched areas in Fig. 0] If 
the x 2 -minimum is bigger than 4 for both types of mass orderings, the corresponding values 
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of |<m>|° s and S are not consistent at 2cr within the assumed uncertainties. Such a 
situation (shown as the shaded regions in Fig. HJ) could either result from systematical ef- 
fects not taken into account in the cosmological data, or can indicate that some mechanism 
beyond the light Majorana neutrino exchange is operating in (/3/3)ov-decay [81]. We also test 
whether the data are consistent with hierarchical spectra, i.e. for each sign of Am^ we test 
whether x 2 { m o = 0) < 4. These regions are below the dashed lines in Fig. 01 within the 
corresponding hatched area. 13 

For experimental errors corresponding to o~pp = 0.03 eV and <te = 0.1 eV adopted in 
the upper row of plots in Fig. IH no distinction between sgn(Am| 1 ) > and sgn(Am| 1 ) < 
(i.e. normal and inverted ordering) is possible. However, some information can be obtained 
on whether the spectrum is hierarchical for a given sg^Am^). In particular, the data 
from cosmology increase the ability to distinguish between IH and QD spectra in the case 
of sgn(Am| 1 ) < if |<m>| ~ 0.1 eV is observed. This situation corresponds to the case 
indicated by the hatched region in Fig. El where (/3/3)oi/-decay alone can only rule out the 
normal hierarchical spectrum. 

For the more demanding experimental precision of app = 0.01 eV and <Te = 0.05 eV, 
used in the lower row of plots in Fig. HI a new possibility to distinguish between normal 
and inverted mass ordering appears. If the {(3(3) 0! ,-decay data give, e.g. a value of |<m>| 
in the interval (0.04 - 0.07) eV for NME uncertainty factor F = 2, and the cosmological 
observations yield an upper bound S obs < 2a s = 0.1 eV, the IH spectrum can be established 
at 2a C.L. This is a qualitatively new method to determine the mass ordering, emerging from 
a synergy between data from (/5/5)ov-decay experiments and cosmology. Using (/?/?) 0l /-decay 
data alone it is possible, in principle, to rule out the inverted ordering if |<m>|° s < 0.01 eV 
and the error opp is sufficiently small. However, as we have shown, the required error 
is exceedingly small: opp ^ few x 10~ 3 eV. In contrast, the conclusion following from a 
combination of (/3/?)oi,-decay and cosmological data is based on i) the observation of a value 
of |<m>| compatible with those predicted for the IH spectrum, and ii) an upper bound 
on mo from cosmological data such that the region where the predictions for | < m > | in the 
cases of normal and inverted hierarchies merge, can be excluded (see Fig. Note that this 
possibility remains even if a factor 3 uncertainty in the NME is taken into account. 

4.5 Constraining the Majorana CPV Phases 

The possibility of establishing CP-violation due to Majorana phases through the ob- 
servation of (/3/3)o;,-decay has been studied previously in [39,41,58]. In the following we 
re-consider this problem by applying the x 2 -method described in Section 14.21 We discuss 
first the case of QD spectrum and later consider the possibility of spectrum with inverted 
hierarchy, sgn(Am^) < 0, and |<m>| obh lying in the IH range of (0.015 — 0.05) eV. As is 
clear from eqs. (fTKJ) and (f2T)|) and the discussion in Section 01 the dependence of | < m > | on 
the phase a^i is suppressed by the small value of sin 2 #13. Therefore, we concentrate on the 
determination of the phase a^i, while the dependence on a^i is taken into account implic- 

13 We adopt the convention to determine the compatibility of the (/3/3)oiy-measurement and cosmological 
data by evaluating the x 2 for 1 degree-of-freedom. A motivation for this convention is provided by the 
so-called parameter-goodness-of-fit method discussed in [82]. In the context of that method the single 
degree-of-freedom corresponds to the one parameter, mo, which is comon to |<m>| and E. 
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Perfectly known NME Factor 1 .5 uncertainty from NME Factor 2 uncertainty from NME 




allowed range for a 21 at la and 2a (1 dof) 
observed \<m>\ = 0.3 eV I I \<m>\ and £ inconsistent at 2a 

sin 2 13 = O± 0.002, sin 2 12 = 0.31 ± 3%, Am^ = 8xl0~ 5 ± 2%, Am^ = 2.2xl0" 3 ± 3% 

Figure 5: Allowed range for the Majorana phase «2i at 1 a C.L. (dark-gray/red regions) 
and 2 a C.L. (medium-gray /blue regions) for |<m>| = 0.3 eV as a function of the observed 
value of E. The shown results are obtained for two sets of assumed errors in the observed 
|<m>| and E (rows of panels) and three values of the NME uncertainty factor F (columns 
of panels). For values of the parameters in the regions between the dotted lines, Majorana 
CP-violation can be established at 2cr. 

itly by minimising the \ 2 with respect to it. For the oscillation parameters we will assume 
throughout this subsection uncertainties at the few percent level (precise numbers are given 
in the figures). Such a precision can be reached in upcoming oscillation experiments. In 
Fig. El we show the allowed range for the Majorana phase c^i for \<m>\ ohs = 0.3 eV as a 
function of the observed mean value of E. The la (2a) range is obtained by the condition 
Ax 2 ( a 2i) = x 2 ( a 2i) — Xmin — 1 (4)- Since the allowed range is determined by A% 2 with 
respect to the x 2 -minimum, there is always an "allowed region", irrespectively of whether 
E obs is consistent with the adopted value of |<m>| obs . We indicate in Fig. El the region 
where the "results" of (/?/3)o^-decay experiment and cosmological observations are inconsis- 
tent (Xmin — 4) by the light shading. Majorana CP-violation can be established if both 
«2i = and «2i = 7r can be excluded. The relevant regions are indicated by the horizon- 
tal dotted lines in Fig. One observes that for app = 0.03 eV and <j£ = 0.1 eV used 
in the upper row of plots, already an uncertainty of a factor of 1.5 in the NME makes it 
practically impossible to establish CPV. Our results show, in agreement with the results 
of the previous studies [41,58], that establishing Majorana CP-violation due to «2i is very 
challenging: the errors in the observed |<m>| and E should not exceed approximately 10% 
and the NME has to be known within a factor F ^ 1.5. Although establishing Majorana 
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uncertainty in \<m>\ from NME 



Y///A data consistent with a = % 
I I \<m>\ and £ inconsistent at 2a 
sin 2 B = O± 0.002, Am^, = 8xl0~ 5 ± 2%, Am^ 

Figure 6: Constraints on the Majorana phase a±2 at 95% CL. from an observed |<m>| obs = 
0.3 eV and (cosmological) "data" on E, as a function of the NME uncertainty factor F. 
Shown are the regions in which i) the "data" are consistent with one of the CP-conserving 
values «i2 = or n (hatched), ii) E obs is inconsistent with |<m>| obs = 0.3 eV (light-shaded), 
and hi) Majorana CP- violation is established (red/dark-shaded). The results are presented 
for three values of sin 2 6 & within the currently allowed range (columns of panels), and for 
two choices of the experimental accuracies for |<m>| and £ (rows of panels). 



1SSSS1 data consistent with a 71 = 
I " 1 CP violation established at 2a 
= 2.2xl0" 3 ±3% observed \<m>\ = 0.3 eV 



CPV would be a very difficult task, it could be possible to exclude a certain fraction of the 
full parameter space of the phase «2i 14 by using the data on |<m>| and E. In particu- 
lar, in many cases it could be possible to exclude one of the CP-conserving values of oc^i, 
a2i = or ct2i = 7r, corresponding to specific relative CP-parities of the neutrinos vi and 
z/ 2 . The sensitivity to «2i depends significantly on the value of the mixing angle 9 & [41,58]. 
As is discussed in Sec. El (see eq. (J22)), for fixed m the allowed range of |<m>| is given 
by m cos 29 & < |<m>| < m . Therefore, the allowed range increases for smaller values 
of cos 28 Q , which makes it easier to exclude the extreme values of |<m>|, corresponding 
to the CP-conserving configurations. This effect is clearly shown in Fig. where the three 
columns of panels correspond to different values of sin 2 ^q. We use the current best fit point 
as well as values between the present 2a and 3a limits. In Fig. |B] we assume an observation 



14 This situation is similar to the case of the determination of the Dirac CP phase S by long-baseline 
oscillation experiments, see [83] for a recent discussion and a list of references. 
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of |<m>|° s = 0.3 eV, which implies a QD spectrum. Adopting representative values for the 
experimental accuracies on |<m>| and E, and scanning values of S obs and the uncertainty 
in the NME, we test first whether S obs and |<m>| obs are consistent at 2a CL. (Xmm — 4)- 
If they are, we test whether the phase ct2i is consistent with the CP-conserving values and 
7i (x 2 (o ; 2i = 0,7r) < 4). These two cases are marked by the hatched regions in Fig. El The 
double-hatched areas, where the regions for ct2i = and tt overlap, correspond to the worst 
situation, where no information on a^i can be obtained and the full range [0, tt] is allowed 
by the data. If the values of E obs and |<m>| obh are consistent and both CP-conserving so- 
lutions for «2i can be excluded, Majorana CP- violation can be established at 95% C.L., and 
we indicate the corresponding regions in red/dark-shading. Clearly, establishing Majorana 
CP-violation becomes possible only under rather specific conditions. For sin 2 ^ 0.31 and 
~ 10% errors in the measured |<m>| obs and S obs (upper middle and right panels in Fig.|Hl), 
the NME has to been known to better than within a factor of 1.5. For smaller values of 
the errors, <jpp ~ 0.01 eV and cte ~ 0.05 eV, Majorana CP-violation could be established 
even for F = 2 (lower middle and right panels). If, however, sin 2 Q = 0.25, the NME 
uncertainty has to be small, F < 1.5 and the indicated high precision in the measurement of 
|<m>| obh and S obs has to be achieved. Finally, the Majorana phase has to have a value 
approximately in the interval ~ (tt/A — 37r/4). 

Consider next the possibility to establish Majorana CP-violation assuming that the v 
mass spectrum is known to be with inverted ordering, sgn(Am\) < 0, and that the observed 
value of |<m>| lies in the IH region of a few x 10~ 2 eV. Knowing that this would require a 
rather precise measurement of | < m >| , we use for the experimental error on | < m> | the value 
a /3/3 = 4 x 10~ 3 eV. For the sum of neutrino masses £ we adopt the error cr^ = 4 x 10~ 2 eV. In 
Fig. [7]we show the sensitivity to Majorana CP- violation for three representative mean values 
of |<m>| from the IH region, |<m>| obs = 0.018; 0.032; 0.047 eV, and two mean values 
of sin , sin 2 = 0.31; 0.38. The allowed regions in all panels of this figure are bounded 
from below by a straight line at S obs = 0.014 eV: below this value £ becomes inconsistent 
with the adopted value of |Am^]- 

The upper row of panels corresponds to the present best fit point of sin 2 . For 
|<m>| obh = 0.018 eV (left panel), the 2a interval of allowed values of |<m>| always 
includes the minimal value of |<m>| for the IH spectrum (see eq. 1)1 Mj) ): |<m>| min = 
a/| Am^| cos 200. The latter corresponds to the CP-conserving value «2i = 7T. Thus, for 
the chosen values of |<m>| obs , app, a-% and sin 2 , it is impossible to establish Majorana 
CP- violation. However, if F < 3, it will be possible to conclude that a^i has a nonzero 
value, Q!2i 7^ 0. Thus, if CP is conserved, the neutrinos i>\ and cannot have the same 
CP-parities. 

The value |<m>|° s = 0.032 eV adopted in the middle panel corresponds to the case 
when |<m>| (with the experimental uncertainty included) satisfies |<m>| min < |<m>| < 
|<m>| max , |<m>| max being maximal value of |<m>| predicted in the case of IH spectrum 
(see eq. (fIH)l ). |<m>| max = y/\ Am\\. If an upper bound on m is provided by a constraint 
on E, Majorana CP- violation can be established, as evident from the red/dark-shaded region 
in the panel. If the observed value of S becomes too large, the "data" becomes consistent 
with the "upturn" of the IH-branch (see Fig. [TJ, which implies that the CP-conserving value 
CU21 = rr is allowed. For S obs < 0.2 eV and uncertainties in the NME F > 1.5 - 2, the "data" 
become consistent with |<m>| max , i.e. with c*2i = 0. 
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sin 2 13 = ± 0.002, Am* = 8xl0~ 5 ± 2%, Am^ = -2.2xl0" 3 ± 3%, a pp = 0.004 eV, o £ = 0.04 eV 



Figure 7: Constraints on the Majorana phase a 12 at 95% CL. for the inverted mass ordering 
from observed values |<m>| obh = 0.018, 0.032, 0.047 eV and (cosmological) "data" on S, 
as a function of the NME uncertainty factor F. Shown are the regions in which i) the 
"data" are consistent with one of the CP-conserving values «i2 = or tt (hatched), ii) E obs 



is inconsistent with |<m>| obs (light-shaded), and iii) Majorana CP- violation is established 
(red/dark-shaded). The upper (lower) row of panels corresponds to sin 2 9 Q = 0.31; (0.38). 



For the third representative value of |<m>|° s = 0.047 eV (right panels), |<m>| max lies 
in the 2a interval of allowed values of |<m>|. This implies that Majorana CP- violation can 
only be established if the "data" on £ constrains tuq precisely in the range corresponding to 
a neutrino mass spectrum with partial inverted hierarchy, such that neither |<m>| max (i.e. 
the horizontal branch at a.21 = 0), nor the upturn of the lower (minimal) branch at «2i = 71 , 
are compatible with the "data" . As can be seen in Fig. [?l this is marginally possible for 
sin 2 9 Q = 0.31 (upper row), but some window exists for sin 2 6* = 0.38 (lower row) if the 
NME uncertainty factor F < 2. 



5 Conclusions 



In the present article we have reanalysed the potential contribution of future (/3/3)ou- 
decay experiments to the studies of neutrino mixing. We have considered 3-v mixing and as- 
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sumed massive Majorana neutrinos and (/3/3) o^-decay generated only by the (V — A) charged 
current weak interaction via the exchange of the three Majorana neutrinos. In this frame- 
work we investigated which information can be obtained from a measurement of the effective 
Majorana mass |<m>| i) on the type of neutrino mass spectrum (NH, IH, QD, etc.) ii) on 
the absolute scale of neutrino masses, and iii) on the Majorana CP-violating phases. As in- 
put in the analysis we used the results of recent studies of the prospective precision that can 
be achieved in the future measurements of neutrino oscillation parameters on which | < m > | 
depends. We performed a \ 2 analysis taking into account experimental and theoretical er- 
rors, as well as the uncertainty implied by the imprecise knowledge of the corresponding 
nuclear matrix element (NME). 

We show how the possibility to discriminate between the NH, IH and QD spectra depends 
on the mean value and the experimental error of |<m>|, and on the NME uncertainty. 
Furthermore, we combine the information on |<m>| from a (f3f3) o v -dec&y experiment, with 
a constraint on the sum of the neutrino masses, S, which can be obtained from cosmo logical 
observations. In this case, we investigate the role of the accuracies on |<m>| and S, as 
well as on the NME uncertainty, in determining the type of neutrino mass spectrum. The 
constraints on Majorana CP-violation phases in the neutrino mixing matrix, that can be 
obtained from a measurement of | < m > | and £ in the cases when i) the observed | < m > | ~ 
few x 1CT 1 eV (QD spectrum), and ii) sgn(Am^) < and the observed |<m>| ~ few x 
1CT 2 eV, are also analyzed in detail. We have estimated the required experimental accuracies 
on both types of measurements, and the required precision in the NME permitting to address 
the issue of Majorana CP-violation in the lepton sector. 

Our results show that, in general, getting quantitative information on the neutrino mass 
and mixing parameters from a measurement of the (/^o^-decay half-life is rather insensitive 
to the errors on the input neutrino oscillation parameters as long as the errors are smaller 
than ~ 10%. However, constraints on the absolute neutrino mass scale, on the type of 
neutrino mass spectrum and on the Majorana CP- violation phase one can obtain depend 
critically on the measured mean value of |<m>| (and £), on the precision reached in the 
measurement of |<m>| (and S), and on the uncertainty in the knowledge of the value of the 
relevant ((3(3)o u -decay nuclear matrix element. The most challenging of these physics goals 
is obtaining quantitative information on Majorana CP-violation phases. The sensitivity 
to the latter depends crucially also on the value of sin 2 # . Establishing Majorana CP- 
violation using data on |<m>| and X in the case of QD spectrum, for instance, would 
require for sin 2 # = 0.31, a ~ 10% (or smaller) errors in the measured |<m>| obs and S obs 
and knowledge of the relevant NME with an uncertainty corresponding to a factor F < 1.5. 
For smaller values of the errors and/or larger values of sin 2 # , say sin 2 # = 0.38, it could 
be possible to obtain evidence of Majorana CP-violation at 2a C.L. even for F = 2. If, 
however, sin 2 6 Q = 0.25, exceedingly high precision in the measurements of |<m>| and E 
and NME uncertainty smaller than 1.5 is required. In all these cases the Majorana phase 
0:21 has to have also a value approximately in the interval ~ (7r/4 — Stt/A). 

Future (f3j3) 0u - decay experiments have a remarkable physics potential. They can establish 
the Majorana nature of neutrinos with definite mass v y If the latter are Majorana parti- 
cles, the (/9/9)o^-decay experiments can provide constraints on the absolute scale of neutrino 
masses and on the type of neutrino mass spectrum. They can also provide unique informa- 
tion on the Majorana CP- violation phases present in the PMNS neutrino mixing matrix. The 
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measurement of |<m>| (and S) with sufficiently small error and sufficiently precise knowl- 
edge of the values of the relevant (/3/3)o w -decay nuclear matrix elements (F < 2) is crucial for 
obtaining significant quantitative information on the neutrino mass and mixing parameters 
from a measurement of (/3/3)oi,-decay half-life. The remarkably challenging physics goal of 
getting evidence for Majorana CP-violation in the lepton sector could possibly be achieved 
only if in addition sin 2 9 Q and the Majorana CP-violation phase ol<i\ have favorable values. 
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